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Abstract 

Exact evaluation of < Tt > is here performed for real symmetric matrices S of 
arbitrary order n, up to some integer p, where the matrix entries are independent 
identically distributed random variables, with an arbitrary probability distribution. 
These expectations are polynomials in the moments of the matrix entries ; they pro- 
vide useful information on the spectral density of the ensemble in the large n limit. 
They also are a straightforward tool to examine a variety of rescalings of the entries 
in the large n hmit. 



1 Introduction 

In the past decades the study and applications of random matrix theory in a variety of 
fields of theoretical physics has witnessed a tremendous expansion. In most cases one 
is interested in universality properties of the model, that is properties of the model, in 
the hmit of infinite order of the matrices, which are independent of the detailed form 
of probability distribution of the random variables. 

Two main classes of models were studied: 

1) Models with independent entries. In this class of models the (independent) ma- 
trix entries are independent identically distributed random variables, that is the joint 
probability density factorizes. 

Several matrix ensembles with real or complex entries were studied. From the sim- 
plest, where no structure of the matrix is considered, to the sparse matrices, see 
for example [2], [3], [5] ; the laplacian matrices (where the diagonal elements are 
Sjj — ~ X^fc-i^n kytj '^i.fe)' ®®® ^'^'^ example [6], [7] ; the band matrices, see for exam- 
ple [1], [12], [13], [15], [14] ; the block band matrices, see for example [17], [18]. 

2) Unitary-invariant models. The three classic models are real symmetric matri- 
ces, complex hermitian matrices or matrices where the entries are real quaternions. 
Here the joint probability density of the random variables is a function invariant un- 
der unitary (real or complex) transformations. The most studied models have the 
joint probability density in the form of a Boltzmann factor which, for real symmetric 
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matrices, has the form 

p(si,i , si,2 , .., s„,„) = e (1.1) 

where V(x) is a polynomial in the variable x. The important case where V(x) = de- 
scribes the gaussian orthogonal ensemble, which is a model also of independent entries. 

Further important matrix ensembles like the Wishart matrices were studied with 
probability density which often do not belong to the above two classes. 

Techniques to study the two classes of models are quite different. In the second 
class, by using the unitarity invarianco, the model is rewritten in terms of the eigen- 
values of the random matrix. Powerful analytic tools are available, like the technique 
of orthogonal polynomials. It was found that while the limiting spectral density is 
sensitive to the details of V{x), several important features related to the local statis- 
tics of eigenvalues exhibit universality properties. 

In the models of the first class few analytical tools are available and often the 
limiting spectral density is obtained by numerical evaluation from samples of large 
matrices. The most famous result is the universality property of the spectral density 
proved by E.Wigner [1] : with suitable assumptions (which will be reviewed later) the 
limiting spectral density is the semi-circle, independent from the probability distribu- 
tion of the entries of the random matrix. 



In this paper few simple ensembles of real symmetric matrices of the first class will 
be analyzed : the {S} ensemble of real symmetric matrices , the closely related {Sq} 
ensemble ( that will be defined in next section), the random sparse matrices and the 
bi-diagonal matrices ( the simplest type of band matrices) . 

Let {S} be the ensemble of real symmetric matrices of order n. Any matrix S G {S} 
is defined by n{n + l)/2 real independent matrix elements , Sjk- The joint probability 
density factorizes 

P , Sl,2 , .., S«,n) = (1-2) 

i<j 

Let X denote any of the independent entries of the random matrix, p(x) its prob- 
ability density, and < x'' > its moments 

£H=<.'>=/^.',(.) fl.3) 

I shall only assume that the moments of any degree exist and in the next section the 
expectations E [TrS''] are exactly evaluated as polynomials in the moments < a;*' > for 
any value of the matrix size n and for p up to 7 for the ensemble {S}, up to p = 8 for the 
related ensemble {5*0}, any p for the bi-diagonal matrices. After a proper n-dependent 
rescaling, the evaluation of E [TrS^], give the moments of the limiting spectral density, 
if E(i) = 0. Symbolic manipulation software is efficient in sorting out words of given 
type and counting them. Then the exact evaluation in the equations (2. 10), (2. 12) 
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may be obtained by such software, by using matrices with symbohc entries of several 
enough orders, say n = 4,5, 6, 7, 8, 9 enough to determine the polynomials in the n 
variable. However computer memory limits a brute force approach, when the number 
of words is a few hundred thousands. In the next section a simple description is given 
in terms of (non-markov) paths which encode all the relevant information. I found 
useful to combine such manual counting with the symbolic software counting. 

The evaluation of the moments of the spectral density, by counting classes of paths 
on the graph associated to the matrix has long history, see for example [1], [8], [9], [3], 
[5]. One may sec [10], [11] for different methods. Usually the aim was to evaluate the 
number of classes of paths for large n and large p, with p. 

The finite p exact evaluation provided in the next section will probably be useful 

for two purposes: 

a) in several models of sparse matrices or laplacian matrices the limiting spectral den- 
sity is determined numerically by extrapolating evaluations for several values of n. In 
such cases, the knowledge of the first few moments of the spectral density, provided 
here, may add useful information. 

b) the results (2. 10), (2. 12) allow a transparent analysis of the n-dependent rcscaling 
preliminary to the n ^ oo limit. Then it becomes evident why universality of the 
limiting spectral density is obtained in the case of Wigner scaling and not in the case 
of sparse matrices scaling. This will be shown in the third section and the possibility 
of different scalings will be considered. 



2 Enumeration of paths 

It is well known that to any square matrix A = {aij}, n x n one may associate a 
graph with n vortices and a directed edge from vertex i to vertex j , corresponding 
to the matrix element Uij. One can the visualize the evaluation of matrix elements of 
powers of j4 as sum of paths on the graph. For example 01,404,202,3 corresponds to the 
three steps path on the graph from vertex 1 to vertex 3 going first to vertex 4 then to 
vertex 2. Similarly the matrix element (^^)i,3 = ^.^j^ ^aijaj,kak,3 corresponds 
to the sum of paths on the graph from vertex 1 to vertex 3 visiting all possible 
intermediate vertices j, k. 

Wc begin by considering the ensemble So of real symmetric matrices of order n, 
with vanishing entries on the diagonal, Sjj = for j = l,2,..,n. Because Sij = Sj,i 
instead of drawing two opposite directed edges between the vertices i and j, one draws 
just one non-directed edge. The graph associated to the generic matrix So is then the 
complete graph with n vertices, shown in the left side of Fig.l 
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Fig.l 

The right side of Fig.l shows the four step path on the graph visiting the se- 
quence of vertices {1, 3, 5, 4, 1}. It corresponds to the product si,3 33,5 55,4 S4,i which 
contributes to Tr 5*0. The same path is depicted on the left side of Fig. 2, where the 
horizontal axes represents time, or number of steps. This type of lattice path is usually 
called an excursion or a bridge respectively if the path is constrained in the upper half 
plane or it is not constrained, see for instance [19]. 



The n{n — l)/2 independent entries of a matrix of the ensemble {^o} are inde- 



pendent identically distributed random variables, 
fe-moment of the probability density. Then 

E [si,3 S3, 5 S5,4 S4,i] =< V >■* 



Let us denote with < 



> the 



(2.1) 



Of course the same contribution is obtained from any other four steps path from vertex 
one to itself, such that the four entries Sij are all different. Because of the symmetry 
of the complete graph, the same is true for any path beginning and ending at a vertex 
different from one. One is led to a more efflcient representation of paths which keeps 
only the relevant information. 





Fig.2 



The paths on the left side of Fig.2 is redrawn on the right side of Fig.2, again with 
time on the horizontal axis. It reads : start from an arbitrary vertex, next a different 
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vertex, next a vertex different from the two previous ones, next a vertex different from 
the three previous ones, finally the first vertex. It is trivial to evaluate the number of 
paths of this kind and obtain the contribution 

n(n - l)(n - 2){n - 3) < v >^ 

I shall call reduced paths the graphs, like the right side of Fig. 2, which correspond to 
a multitude of paths where the visited sites are listed, like in the left side of Fig.2. 
To evaluate < Tr Sq >, only four reduced paths , shown in Fig. 3, are needed. Their 
contribution is 

<Tr5'o> = n{n - 1) < V* > +2n{n - l){n - 2) < >^ + 

+ n(n - - 2)(n - 3) < >■* (2.2) 




Fig.3 

From now on I shall only refer to these type of paths. 

I shall call closed paths those where the final vertex is the same as the beginning 
vertex. Open paths have a final vertex different from the beginning vertex. Only 
closed paths contribute to the evaluation of Tr Sq, but open paths are equally impor- 
tant because all the closed paths with p steps are obtained by considering all the open 
paths with p — 1 steps and adding one further edge connecting the last vertex to the 
first vertex (closing the path). 

Let Sc{n,q) denote the number of closed paths with n steps and such that the 
maximal altitude reached during the path is q. It means that the path visited a total 
of q distinct vertices of the complete graph, including the initial vertex. Similarly 
Soin, q) will denote the number of open paths with n steps and such that the maximal 
altitude reached during the path is q. The previous remark translates into 
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Sc{n,q) = So{n - l,q) 



(2.3) 



One may evaluate the total number of paths , that is the open plus closed paths 
for any number of steps. 

Consider a path with n — 1 steps and denote byj/jt ,fc = 0,l,---,n— 1, the altitude 

of the path at time k. Let q— maximum{j/fc}. One further step may be added to the 
path, by connecting yn-i with where j/n is any of the q+1 integers 0, 1, 2, .., g + 1, 
excluding j/„ = j/„_i. Of course there is no path where n < q. 

Let S{n,q) = Sc{n,q) + 5'o(n,q) be the total number of paths with n steps and 
maximal altitude q. The previous construction of paths of n steps from paths with 
n — 1 steps , translates into the equations 



S{n,q) = qS{n-l,q) + S{n-l,q-l) 

S{n, q)=0 if n< q ; S{n, 1) = 1 (2.4) 

The positive integer numbers S{n, q) arc known as Stirling number of the second 
kind. Their sum over the values of q are known as the Bell numbers B(n), see for 
example [21], [20]. Useful formulas are 

S{n,q) = \,Y.i-^y-{l)r^ > 

n 

B{n) = ^5'(n,g) , 

q=l 

n 

B{n+1) = ' 

r=0 
oo 

exp(e^-l) = ^^x"B(n) (2.5) 

n=0 

For example the total number of paths with four steps is B{4) = 15 and 5(4, 1) = 1, 
5(4, 2) = 7, 5(4,3) =6, 5(4,4) = 1. 

The first few Bell numbers are 1 , 2 , 5 , 15 , 52 , 203 , 877 , 4140 , ■ • • 

Because of eq.(2.3), 5c(n) = 5c,(n — 1) = B{n — 1) — 5c(n — 1), then the first few 

numbers of closed paths are 

5c(2) = l , 5e(3) = l , 5c(4) = 4 , 5c(5) = 11 , 5c (6) = 41 
5c(7) = 162 , 5c(8) = 715 (2.6) 

By substituting all the non- vanishing entries of the matrix So with one, a circulant 
matrix is obtained where it is easy to evaluate 

TV5o^^ = (n-l)[(n-l)^^-^ + l] , 
Tr5o''+' = (n-l)[(n- 1)^^-1] (2.7) 
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One obtains a new and easy determination of the number of closed paths Sc{n, q) by 
comparing the above numbers with the contribution of the graphs. For instance, in 
the case of paths with 8 steps, by solving the equation 

8 

(n - 1) ((n - I)'' + l) = ^ n{n - l)(n - 2) • • • (n - g) Sc{8, q) (2.8) 

9=1 

We obtain 

Sa(8,l)^l , S'c(8,2)=42 , S'c(8,3)=231 

5c(8,4)= 294 , S'c(8,5) = 126 , S'c(8,6) = 20 

8 

5c(8,7) = l , 5^Sc(8,9) = 715 (2.9) 

9=1 

Finally one draws the graphs and reads their contributions. It is a time consuming 
activity, helped by symbolic manipulation software. I report here the evaluations up 
to E [Tr^o], which results from 715 graphs. 



- < TrSo > 
n 

- < Tr5o > 
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< Tr So > = 
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l)<v^> 

l){n~2) <v 

!)<«''> +2(n - l)(n - 2) < >^ + 
l)(n-2)(n-3) <w >" 

- l)(n- 2) < >< V >^ + 

- l)(n - 2)(n - 3) < >< w >^ + 
l)(n-2)(n-3)(n-4) < v >^ 

1) <v'^ > +6(n - l)(n - 2) < V* >< > + 

l)(n-2)(5n- 11) < >^ + 

- l)(n - 2)(n - 3) < >< t; >^ + 

- l)(n - 2)(n - 3)(2n - 5) < >< t; + 
l)(n - 2f{n - 3)(n - 4) < w 



i<TrS'^> = 7(n- l)(n-2) < w'^ >< w >^ +14(n- l)(n-2) < t;^ >^< v > + 



l)(n - 2)(n - 3) < w >< w >^ + 
- l)(n - 2)(n - 3) < i;^ >< >< v >^ + 
l){n - 2){n - 3)^(3n - 8) < >^< v >^ + 



+ 7{n - l){n - 2){n - 3){n - 4) < >< V + 
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+ 7(n - l)(n - 2f{n - 3)(n - 4) < >< v >^ + 
+ (n - l)(n - 2)(n - 3)(n - 4)(n^ - 4n + 2) < t; 

-<TrS^> = (n- 1) < > +(n- l)(n-2)[8 < >< > +6 < + 
n 

+ {n - l){n - 2f28 < >< v'^ >^ + 

+ (n - l)(n - 2)(n - 3)[8 < ■y'^ >< « >^ +20 < >^< v >^] + 

+ (n - l)(n - 2f{n - 3)48 < v'^ >< >< v + 

+ 4(n - l)(n - 2)(n - 3){2n - 3) < >< w >^ + 

+ (n- l)(n- 2)(n- 3)(14n- 19) < -y^ >^ + 

+ 2{n - l)(n - 2)(n - 3)(14n^ - 66n + 51) < >^< v >" + 

+ 8n(n - l)(n - 2)(n - 3)(n - 4) < >< w + 

+ 8(n - l)(n - 2){n - 3)(n - 4)(n^ - 3n - 2) < >< v + 

+ (n- l)(n - 2)(n - 3)(n - 4)(n - 5)(n^ - n - 4) < w (2.10) 

M.Bauer and O.Golinelli, in their analysis of the incidence matrix for random 
graphs [4] , evaluated traces of powers of a random matrix by evaluating normaUzed 
fc-plets. These arc sequences of k positive integers {vi,V2, • • • , Vk) such that : 

1) Vl ^ V2 , V2 7^ V-i, ■ ■ ■,Vk-l 7^ Vl 

2) if > 1 is in a normalized /s-plets, then also Vpi = U/j — 1 is in the same normalized 
fc-plet and /3' < /3. 

The normalized /c-plets are in one-to-one correspondence with the reduced closed paths 

described in this paper : the sequence of the k integers of a normahzcd fc-plct is the 
sequence of the y-coordinatcs {yo = 1, y2, ■ ■ ■ , Vk-i}, of a reduced closed paths with 
k steps. For example, the four reduced paths in Fig. 3 correspond to the four nor- 
malized 4-plets {1,2,1,2} , {1,2,1,3} , {1,2,3,1} , {1,2,3,4}. 1 prefer the graph 
representation because 1 find it easier to read the repeated steps and it allows several 
generalizations for random matrices with diagonal elements not equal to zero, which 
may have the same probability distribution of the other entries, as in the next section, 
or a different probability distribution, as it occurs in laplacian matrices, not described 
in this paper. 



From the ensemble {5*0} to the ensemble {S}. 

The graph corresponding to a matrix of the ensemble {S}, of order n, is the 
complete graph with n vertices, with a loop added to each vertex. It is depicted 
in Fig. 4. 
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Fig.4 



The evaluation of E [Tr S^] may be performed by considering the closed 
graphs with p steps where any number of horizontal steps are now allowed at 
any vertex. A simple example will show how the new reduced paths may be 
obtained from the reduced paths contributing to the evaluation of E [Tr Sq] , 
fc = 0, 1, • • • ,p , (where horizontal steps are missing) by adding the p ^ k hori- 
zontal steps in all possible ways. 

For example to evaluate E [TrS"*] we consider the two paths, in Fig. 5, respec- 
tively of 2 and 3 steps. 




Let us first consider the graph of 3 steps. It has 4 vertices and we obtain 4 
closed graphs contributing to E [Tr S'^] by adding one horizontal step in four 
possible ways , as it is shown in Fig. 6. Each of the 4 graphs corresponds to the 
contribution 

n{n - l)(n - 2) < v 
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Fig.6 

Next one adds two horizontal steps to the graph of 2 steps in Fig. 5, in 
6 possible ways, thus obtaing the 6 graphs shown in Fig. 7. In the first 4 of 
them , the two horizontal steps are at the same vertex and their contribution is 
4n(n — 1) < >^. The last 2 graphs have the two horizontal steps in different 
vertices, then their contribution is 2n(n — 1) < v'^ >< v >^ 





Fig. 7 

Finally one adds the contribution of one graph with four horizontal steps , that 
is n < u'* >. By collecting the new terms one obtains 

E [Tr S"*] = n<v'^ > +4n{n - 1) < >2 +2n{n - I) < v"^ >< v y"^ + 
+ 4n(n - l)(n - 2) < t; +E [Tr S^] = 

= < > +2n{n - I) < v'^ >< v >^ +2n^(n - 1) < i;^ >^ + 
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+ n{n^ - l){n - 2) < V >^ (2.11) 

As a trivial check, by replacing all the moments < >^ 1 , one finds 

E[TtS^] =n*. 

One may then evaluate 











- < Tr5 


> 




<v> 


n 








-<TrS^ 


> 




n Kv^ 


n 








- < TrS^ 


> 




<v^ > 


n 








- <TV5^ 


> 




n <v^ 


n 












+ 


2n{n — 



I) Kv^ ><v > +(n - l)(n - 2) < t; >' 



,2 >2 + 
2 i\/„ o\ ^ .. ^ 4 



+ {n^ - l){n - 2) < V >' 

-<TrS'^> = <v^ > +bin-l) Kv"^ ><v> +bin-l) <v^ Xv^ > + 
n 

+ 5(n - l)(2n - 3) < i)^ >^< i; > +5(n - l)(n - 2) < t;^ >< t; >^ + 

+ .5(n - l)(r?, - 2f < >< v >^ + 

+ {n - l){n - 2){n^ - 2n + 2) < V >^ , 

-<TrS'''> = n <v^ > +6(n- 1) <v^ XV +6(n-l) <v^ Xv'^ XV > + 
n 

+ 3(n- l)(2n+ 1) < >< > + 

+ 15(n- l)(n- 2) < >2< w >2 + 

+ (n - l)(5n^ - 6n - 5) < >^ + 

+ 6(n- l)(n- 2)(n + 3) < -u^ >< i; >^ + 

+ 3(n - l)(n - 2)(2n^ + n - 17) < >< v >^ + 

+ (n - l)(n - 2)(n^ - Sn^ -7n + 23) <v >^ 

-<TrS'^> = <v'' > +7{n-l) <v^ XV > +14:(n-l) <v^ ><v^ > + 
n 

+ 7{n-l) <v^ Xv"^ > + 

+ U{n- l){3n- 4:) <v^ Xv'^ XV > + 

+ 7(n- l)(n-2) < >< i; >^ + 

+ 7{n-l){3n-5) <v^ ><v^ >^ + 

+ 14(n - l)(n - 2) < >^< w > + 

+ 7{n-l){n-2f <v'^ ><v>^ + 

+ 7{n-l){n-2){5n-8) <v^ ><v^ XV >^ + 

+ 35{n-lf{n-2)<v'^>^<v> + 
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+ 7{n-l)in-2){n^ + 1) <v^ XV >* + 

+ 7(n-l)2(n-2)(3n-8) <t;^ >^<t;>^ + 

+ 7{n - l)(n - 2)(n^ - 2n^ - 2n - 2) < v"^ >< v >^ + 

+ (n - l)(n - 2)(n - 3)(n^ -n^ -Wn-l)<v >^ 

(2.12) 

The bi-diagonal symmetric matrices. 



In a lattice in one dimension with n sites and periodic boundary conditions, 
each site has two next neighbours, the lattice has n edges. To each edge {k, k+1) 
we associate the random variable ak,k+i = afe+i,fe- 



We consider the real symmetric bi-diagonal random matrix A 



( 



A = 






Ol,2 











a\,n 


ai,2 





02,3 














02,3 





a3,4 



































^n — l.n 






(2.13) 



A walk on the lattice with p steps may return to the initial site only if p is even. 
The number of such paths ^ is " This is the number of terms contributing 

to i < TvAP >. 

The left side of Fig. 8 shows the periodic one dimensional lattice with n vertices, 
associated to the matrix A in eq.(2.13). The right side of Fig. 8 shows a path 
from site j returning to site j after 14 steps, then contributing to Tr A^^. 
Any path may be seen as composed of paths staying (strictly or weakly) in the 
positive or in the negative half-plane. Paths which do not enter the negative 
half-plane (these are the weakly positive paths), such as the one depicted in 
the left side of Fig. 9 , are known as Dyck paths. Dyck paths have a one-to-one 
correspondence with rooted trees, as indicated in the right side of Fig. 9. It is 
the possible to enumerate the number of steps (up or down) between level k and 
level k+1 for the class od Dyck paths of any length. Next by composing wciakly 
positive and weakly negative Dyck paths, one obtains enumeration of steps for 
the generic paths, see for instance [22]. 

^Because of the periodic boundary conditions, there exist walks on the lattice which return 
to the initial site from the opposite side of the first step, if p > n. Since we are interested in 
the limit n — » oo with fixed p, those walks are here neglected. 



12 



1 2 




Fig.8 



One easily evaluates 





> = 


, 


- < TiA" 1 
n 




n 


> = 


2<v^ > 






- < TiA^ 
n 


> = 


2 < w'' > H 


-4 < >2 




- < TrA^ 
n 


> = 


2<v^ >^ 


-12 < v'^ >< > ^ 




- < Tr A^ 
n 


> = 


2 < > H 


-16 < -y*^ >< -y^ > 4 


12 < -y'' >2 


+ 


32 < >< >2 +8 < >^ 




- <TrAi° 
n 


> = 


2<v'"> 


+20 < >< > - 


f40 <v^ >< 



+ 50 <v^ >< >2 +70 < >2< > +60 < >< >^ + 

+ 10 < >^ (2.14) 




Fig.9 
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Fig. 9 shows the one-to-one correspondence between Dyck paths and rooted 
trees : given the tree, one may follow its external contour, starting at the left 
side of the root. As one moves up along the first branch, one generates two up 
steps of the Dyck path , then moving down along the same branch generates the 
two steps down ; next branch is shorter corresponding to one step up and one 
step down for the Dyck path. Finally, after touring next branch, one reaches 
the right side of the root. 



3 Scalings and/or Universality 

To study the existence of a limiting spectral density in the n — > oo limit, some 

n-dcpendent scaling for the entries of the matrices is needed. 
One may recall that the existence of finite large-n limit of a finite number of 
rescaled traces of powers of the random matrices is neither necessary nor suf- 
ficient condition for the existence of a limiting spectral density. Indeed in the 
case of gaussian matrices when < v > is not equal to zero and Wigner scaling 
(which will be recalled here below) the large-n spectral density of the ensemble 
contains the classical semicircle separated by a small peak, which disappears at 
n = DO, yet produces the large-n behaviour of the traces.'^. A more complex dis- 
tribution with possibly several peaks occurs with the sparse matrices [4], [8], [7]. 

Conversely, a smooth limiting spectral density for the matrix ensemble may 
exist , with unlimited support, and decreasing as an inverse power for large 
eigenvalues. Then only a finite number of moments of the spectral density will 
exist. 

Then the aim of this section is not to derive the limiting spectral density 
from the table of values (2.10), (2.12) for < TrS'P > and < Tr^^ > after proper 
rescaling, but merely to review the effects of two well known scalings and to 
indicate the existence of intermediate scalings which appear to yield the same 
results of the Wigner scaling for every finite p. 

Wigner scaling 

The Wigner scaling consists of two assumptions : 

a) the i.i.d. entries have vanishing expectation value, that is E {sij) =< v >= 
for any j = 1, 2, • • • , n. 

b) the i.i.d. entries are written 

Sij = , ^ij independent on n 

v n 

suppose this is well known to experts, but I have not found a proper reference, then I 
offer the short discussion in the Appendix, by using the addition law for fixed plus random 
ensembles. 
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then E [(sij)'^] =< v'^ >= —^Ck , Ck independent on n 

By inserting these assumptions in the tables (2.10), (2.12) and by keeping only 
the leading order contribution in the limit n — > oo , both tables collapse into 
the very simple result 



< Tr 5^ > 



- < TrS"^ > 
n 

- < TvS^ > 
n 

- < Tr5^ > 
n 



< r > 



2 < >^ +0 - 
n 



5 < ^2 >3 +0 ( - 

n 



14 <^'^>'^+0[- 

n 



(3.2) 



Which are easily seen as the lowest moments of the semi-circle density 
n \n 



(3.3) 



The variance a =< ^ > of the probability density of the random variables 
is the only parameter entering in the limiting spectral density 



G{z) = lim -E 

n— >oo 77, 



Tr- 



z-S 



oo 

= lim -y"_-E[Tr5'f] 



p=0 



P(A) 



fe=0 •' ^ 

TT ^ \ if |A|>2a 



Random graph scaling 
The simplest model of random graph is the ensemble {iSq}, with the probability 
density 



1 , with probability 1/n , , , 

, with probability 1— 1/n , \ - ) 
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Then E [(sij)*^] =< v'' >= 1/n for every k. By inserting this in the tables 
(2.10), (2.12) and by keeping only the leading order contribution in the limit 
n ^ oo , one finds that all the odd moments and therefore the odd powers of 
traces are negligible. The limiting form of the tables become 



lim - < TrS'^ > = 1 
lim - < TrS^ > = 3 

n— »oo fi 

lim - < TrS''^ > = 12 

n^oo Tl 

lim - < TrS^ > = 57 

n— >cx) n 



(3.6) 



These numbers reproduce the asymptotic behaviour in the limit n ^ oo of the 
moments E [iTrS^^] = nip , evaluated by A.Khorunzhy and V.Vengerovsky 
[3] . They also found a recursive relation that evaluates the asymptotic moments 
nip and studied their behaviour at increasing values of p. The same recursive 
relation was obtained and analyzed by M.Bauer and O.Golinelli [4]. 

A very similar scaling is used in models of speirse matrices : 

< {s.j)'' >=< v'' >= ^Cfc (3.7) 

It is easy to see from the tables (2.10), (2.12) that with this rcscaling same 
limit holds for matrices of the two ensembles {S} and {5*0}) and that the odd 
moments vanish 

lim - < TrS'^f+i >= 

n— >cso 77, 

C2 

C4 + 2(C2)2 

C6 + 6C4C2 + 5(C2)^ , 

C8 + 8C6C2 + 6(C4)^ + 28C4(C2)^ + 14(C2)^ , (3.8) 

Intermediate scalings 

In both scalings above mentioned, the Wigncr scaling and the sparse matrix 
scaling, the second moment of the entries is scaled with the same power 

< {Sijf >=< >= -C2 
n 



lim — 


< TrS'^ 


> = 








lim — 


< Tr5^ 


> = 


n — *oo fi 






lim — 


< TrS"5 


> = 








lim — 


< TrS'^ 


> = 


n— >cx) n 
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Every moment < v'' > with fc > 2 , is scaled by the sparse matrix scaling (3.7) 
in a way sufficient to obtain finite asymptotic limits lim„^oo ^ < Tr S'p >. Fur- 
thermore the same scaling suppresses asymptotically every contribution of the 
first moment < v >— ci/n. 

If the first moment < w > is not zero, and the matrix entries have Wigner 
scaling, < v'' >= Ck/n^^"^, see eq.(3.1), several rescaled E[TrS'P] or E[Tr5o] 
diverge in the n — *■ cx3 limit. 1 considered a simple gaussian example in the Ap- 
pendix. Most authors assume the first moment < w >= together with Wigner 
scaling, as I did at the beginning of this section. 

Wigner scaling suppresses asymptotically every contribution < > with 
k > 2 and in this way obtains universality. 

Let us parametrize possible scalings of the moments of the probability dis- 
tribution of the entries 

< is,,, f >=< v"^ -^Ck (3.9) 

Fig. 10 plots the rescaling power a(jj) on the vertical axes. The horizontal line 
a(j}) = 1 corresponds to the sparse matrix scaling ; the line a{p) = p/2, with 
p> 2 corresponds to Wigner scaling. 




Fig.lO 



Any line in the sector between the two lines, that is any scaling with a(p) = 
ap + 1 — 2a with < a < 1/2 and p > 2, supplemented by < w >= 0, would 
obtain finite asymptotic evaluations for ^ < TrS'^' >. Our tables, up to p = 8 
suggest that such asymptotic evaluations would be identical with Wigner scaling, 
provided the first moment vanishes, < v 0. 

A simple example of distribution probability corresponding to intermediate 
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scaling is 

, with probability 1 — ^i_2„ , 
1 



with probability in^--ia > (3.10) 



where < a < 1/2 ; then one evaluates 

< is.,^ >- ( ? ' l ^'^^ ' (3.11) 

I „a(fc-2)+i , if A; even ^ ^ 

4 Appendix 

The effects of a non- vanishing expectation value < v > may be seen in the simple 
model of the ensemble {S} of real symmetric matrices with normal probability 

distribution for the entries. 

The joint probability density factorizes 

P (Sl,2 , Si,3 , .., Sn-l,n) = T[p {Si,j) , P {Sij) = -J—e'^''^-'"^'/'^''^ (4.1) 

The first moments < > are well known 

< Sij >=V , < {Sijf >= CT^ + 

< (sij)^ >= 3a'^v + , < {sij f >= Sa"^ + 6a^v^ + 

< {sijY >= IOSctS + 105crV + 21a'^v^ + 

By inserting these moments, after Wigner rescaling Sij Sij / \fn , in the table 
(2.12), one finds the large n behavior of < TrS" > for p = 1, 2, .., 7. 

-<Tr5> = v/^/n 
n 

-<TrS^> = {a^+v^) , 
n 

-<TrS^> = ni/2^;3 + 3n-i/Vw , 
n 

-<TiS^> = nv^ + 2{a'' + 2v^)a'' + -{a^+Av'')cr^ + o(\ 

n n yw^ 

-<TrS'^> = n' 
n 



3/2^5 ^ 5^1/2^2^3 + ^(^2 + V^^V + O ( 



n \n 

-<TrS'> = n5/2^;7 + 7n3/Vt;^ + 7nV2(3a2+4^;2)aV + 0(l) (4.3) 
n 



18 



The large-n resolvent G{z) of this ensemble may be investigated with the 
addition law of random matrices. [23], [24], [25], [26]. 



Any real symmetric matrix of the ensemble {S} , where < v > is not equal 
to zero, is the sum of a random matrix of the same ensemble, with < v >= 
and the fixed matrix u J, where the J matrix has all entries equal to one. After 
Wigner scaling all the entries of the random and fixed matrices have the further 
scaling factor n~^^'^. 

The resolvent G{z) and its inverse function, for the fixed matrix vJn~^^'^ 

are 

n- 1 1 1 1 



G{z 



n z 



1 + ^vG + v^(l - ^vGf + 4n-i/2^G 

The large— n limit for the resolvent and its inverse function, for the ensemble 
of gaussian matrices with <v>—Q correspond to the Wigner semi-circle 



G{,z) = ±^(z-^/z^-Aa^) 



2<t2 

z = l+a^G (4.5) 

According to the law of addition of random matrices, the resolvent and its 
inverse function, for the ensemble {Sn~^/^}, for large-n have the Pastur form 



G{z) = -V — 



n-l 1 1 

+ -- 



n ^ z — Cj — a'^G n z — a'^G n z — V\/n — a'^G 
3 ■' 



z = ,2^^^ + V"rG+^il-VnrGr + ,n-^n,G^ ^^^^ 

On the right side of eq.(4.6) for the resolvent G{z) , the first term dominates the 
second one for every z finite and n large, whereas the second term is dominant 
for z very close to v^/n. Then a simple approximate form follows by considering 
only one term 



Giz) 



2^ (^z — \/ z^ — 4(1 — l/n)cr2^ z finite and n large 

^(^z- vy^ - - - 4a2/n) z ~ v^/n + ^ 



2a- 

corresponding to the spectral densities 
1 

27ra 



pW ~ 7r-^V4(l-l/n)cT2-A2 , if |A| < 2aV(l- 1/n) , n large 

ZTTCr 
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Here the upper line reproduces the usual semi-circle spectral density, while the 
lower line is a very narrow peak at A = Vsfn. 

Finally the eq.(4.6) for z = z{G) is easily expanded for small values of G, 
the expansion is inverted obtaining the large z expansion of the resolvent, as it 
follows from the addition law 

+ (n^/\^ + 5nV\3^2 ^ 5^-1/2^^2(^2 ^ 2^2)) 1 + 

+ ( + Qna\* + (5^4 ^ ^5^2^2 ^ 9^4)^2 ^ 15^4^2 
V n 

+ (n'/W + 7n^/^a\^ + 7n^^\2v^ + 3a^)a\^) ^ + • • • (4.7) 

By comparing this expansion with the direct evaluation of the large n behaviour 
of E(TrS'^') given in (4.3), we see that the first two terms of the expansion are 
correctly reproduced for every p = 1, 2, • • • , 7. Analogous analysis and conclu- 
sions are abtained for the {So} ensemble. 
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Figure Captions. 

Fig.l. The left side shows the complete graph with N vertices, 
the right side shows a 4 steps path on it. 

Fig. 2. In both graphs of the figure, time is the horizontal axis. 
The left side shows the history of a path on the graph by exhibiting 
the visited sites. The right side is the graph of a reduced path cor- 
risponding to the graph on the left side. 

Fig. 3. The four reduced paths which yield < Tr5'o >. 

Fig. 4. The complete graph with N vertices corresponding to the 
symmetric matrix S, with non-zero entries on the diagonal. 

Fig. 5. Two reduced paths of two and three steps. 

Fig. 6. One horizontal step has been added in 4 possible ways to 
the 3 steps graph of Fig. 5. 

Fig. 7. The six graphs arc generated by adding 2 horizontal steps 
to the 2 steps graph of Fig. 5. 

Fig. 8. The left side represents the graph corresponding to the 
bi-diagonal matrix in eq.(2.13). The right side represents a 14 steps 
paths on the graph beginning and ending at site j. Time is the 
horizontal axis. 

Fig. 9. The one-to-one correspondence between Dyck paths and 
rooted trees. 

Fig. 10. The re-scaling power a{p). The line a = 1 is the re-scaling 
used in random graphs, the line a{p) = p/2 is the Wigner scaling. 
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